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Abstract. A periodic graph in dimension d is a directed graph with a free ac- 
tion of Z d with only finitely many orbits. It can conveniently be represented in 
terms of an associated finite graph with weights in Z d . Here we use the weight 
sums along cycles in this associated graph to construct a certain polytope 
in K. d , which we regard as a geometrical invariant associated to the periodic 
graph. It is the unit ball of a norm on R d describing the large-scale geometry 
of the graph. It has a physical interpretation as the set of attainable velocities 
of a particle on the graph which can hop along one edge per timestep. Since 
a polytope necessarily has distinguished directions, there is no periodic graph 
for which this velocity set is isotropic. In the context of classical physics, this 
can be viewed as a no-go theorem for the emergence of an isotropic space from 
a discrete structure. 



1. Introduction 

Periodic graphs are abstractions of the atomic structure of crystals. A crystal, 
by definition, is a material whose structure consists of a finite-size pattern which re- 
peats periodically in all spatial directions. Taking the crystal atoms as the vertices 
of a graph and the chemical bonds as its edges, one obtains a graph which repeats 
periodically in all spatial directions: a periodic graph. This graph represents the 
chemical structure of the crystal. Therefore, the problem of classifying and enu- 
merating all periodic graphs in three dimensions is of fundamental importance for 
crystallography [2, 8, 9, 19]. Periodic graphs have also been studied in operations 
research [16], spectral graph theory [5], and computer science [4]. Certain gener- 
alizations of periodic graphs also appear in topological graph theory [13, 14, 20]. 

Besides their natural appearance in all these fields, periodic graphs have re- 
cently also been suggested as candidates for the fundamental microscopic building- 
blocks of space [6, 7). In order to get some idea of how physics on a periodic graph 
can look like, let us consider a classical point particle moving on a periodic graph 
r C M. d as follows: the particle moves along the vertices of T in discrete timesteps 
by hopping along one edge per timestep. More precisely, we define a trajectory 
to be a sequence (/ n )neN of vertices f n S T such that for each time n £ N, the 
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positions f n and f n +i are adjacent in T. A vector u e W l is then called a velocity 
vector of T if there is a trajectory / such that 



n— >oo 

Intuitively, this equation means that the trajectory's apparent velocity on the 
macroscopic scale is given by u. It is the trajectory's velocity as seen by a macro- 
scopic observer who is not aware of the fundamental discreteness of T and perceives 
space as a continuum M. d . Note that (1) only makes sense when the limit on the 
right-hand side exists, which can be interpreted as requiring the trajectory to have 
a well-defined constant macroscopic velocity. 

Now a natural question is: given the periodic graph V e M d , what is the set of its 
velocity vectors? In particular, can this set be a Euclidean ball, thereby making the 
macroscopic observer perceive an isotropic space, such that the achievable absolute 
values ||u|| do not depend on the direction u/||u||? This would be a very desirable 
property for the kind of models discussed e.g. in [6, 7]. 

Using concepts from the theory of periodic graphs, we will prove in Theorem 16 
that the set of velocity vectors of any suitably connected periodic graph T is a convex 
polytope in M. d . In particular, it never is a Euclidean ball, and the set of achievable 
velocity vectors cannot be isotropic; see Theorem 25. This is a no-go theorem for 
the emergence of an isotropic space from a discrete structure within the context of 
classical physics. 

From the mathematical point of view, our velocity polytopes are new invari- 
ants of periodic graphs. As witnessed by Proposition 19, they encode the periodic 
graph's large-scale geometry. The velocity polytope as an invariant can be applied 
for example as in Corollary 24, which is a criterion for proving the non-existence of 
translation-invariant maps between periodic graphs ( "morphisms" ) . 



2. Preliminaries 

In this section, we collect some definitions and simple observations. Although 
most of the relevant literature is concerned with the case of undirected graphs [2, 
3, 4], we work with directed graphs, which is more general and has turned out to 
be technically more convenient. 

If A is a finite set, we write \A\ for its cardinality. 

2.1. Graphs and paths. For us, a graph is a directed graph which may have 
loops and multiple edges. A graph G is specified by a vertex set Vq, an edge set 
Eq, a source function sq ■ Eg — > Vq, and a target function to ■ Eq — > Vq- We 
refrain from identifying an edge e with the vertex pair (sg(s), tcie)), since typically 
there will be several edges between sc(e) and tc(e). When the graph G is clear 
from the context, we frequently omit the subscripts and simply write s,t : E — > V 
for the source and target maps in order to avoid unnecessary cluttering. 

A path p in G is a finite sequence of edges p = e± . . . e„, S E, such that 
i(ej) = s(e; + i) for all i = 1, . . . , n — 1. The length \p\ = n of p is its number of 
edges. The empty path is the unique path of length 0. 

A closed path is a non-empty path p = e\ ...e n such that s(ei) — t(e„). If 
no other additional vertex repetitions occur, then p is also called a cycle. By this 
definition, a cycle of length n traverses n distinct vertices. 
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If p = ei . . . e n and p' = e[ . . . e' m are paths such that t{e n ) = s(e[), then p and 
p' can be composed to a new path 

I dof / / 

pp = ei . . . e n e 1 . . . e m . 

It is clear that \pp'\ = \p\ + \p'\. 

A directed graph is said to be strongly connected if there is a path from v to w 
for any pair of vertices (v, w). 

Lemma 1. Let G be a finite directed graph. 

(a) A path p = e\ . . .e n in G of length n > \Vg\ contains a cycle: there are 
indices k and I such that c = e^ek+i ■ ■ ■ s-i is a cycle. 

(b) There are only a finite number of cycles in G. 

Proof. (a) By the pigeonhole principle, there have to be indices I > k 
with t{ei) = s(efc). Choosing I minimal with this property guarantees 
c = ekCk+i ■ • • e; to be a cycle. 

(b) Since a cycle is defined as not having any vertex repetitions besides the 
coincidence between the initial and the final vertex, a cycle can have length 
at most | Vq I . The conclusion follows since there are only a finite number 
of paths of length at most |Vg|- 

□ 

2.2. Definition of periodic graphs. We now turn to the formal definition 
of periodic graphs before discussing their representation by finite weighted graphs. 
For examples, we refer to figure 1 for a basic example in two dimensions and to [2] 
for abundant visualizations of three-dimensional periodic graphs within the context 
of crsytallography. 

Definition 2 (periodic graph). A d- dimensional periodic graph is a graph T equipped 
with a free action of the free abelian group Z d on T , such that T has only finitely 
many *Z d -orbits of vertices as well as edges. 

Let us disentangle what this definition means. First of all, the graph T comes 
with an action of the group 7L d . In additive notation, this means that there are 
given maps 

VrxZ d ^V r , (v,x)^v + x, 
E r xZ d ^E r , (e, x) >->• e + x . 

We think of the vertex v + x as the vertex v translated by the vector x £ Z d , and 
similarly for e + x. In order for these maps to form a Z d -action on T, they need to 
satisfy the group action axioms 

(v + x) + y = v + (x + y) \/v G Vr, x, y G Z , v + = v \/v G Vr , 

(2) 

(e + x) + y = e + (x + y) Ve G E T , x, y G Z d , e + = e Ve G E? , 

as well as be compatible with each other in the sense that source and target of the 
translate of an edge are precisely the translates of the source and target of the edge, 
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sr(e + x) = sr(e) + x , 

(3) 

t r (e + x) =t r (e)+x. 

Furthermore, the action of Z d on T should be free, 

v + x = v + x=$-x = x, e + x — e + x => x = x . (4) 

Finally, there should be only finitely many orbits in Vr as well as in Er under the 
Z d -action. This is an abstraction of the crystallographic property that a unit cell 
of a crystal contains only finitely atoms and chemical bonds. 

Remark 3. Definition 2 is of an abstract combinatorial nature in the sense that no 
embedding of T into M. d is required. However, suppose that r C M. d is a Euclidean 
periodic graph which is translation-invariant in d linearly independent directions. 
After applying an appropriate affine transformation to T, the unit cell of T can 
be taken to be the unit cube [0, l] d , which means that T is translation-invariant 
under the group of integer translations Z d C K d . With Z d acting on T by these 
translations, T is a d-dimcnsional periodic graph in the sense of definition 2. 

There are other uses of the term "periodic graph" in the mathematical litera- 
ture which are not related to the one used here. For example, a plot of a periodic 
functions is a "periodic graph" in a completely different sense. For another in- 
teresting notion of periodic graph which is in no way related to the present one, 
see [11]. 

2.3. Displacement graphs of periodic graphs. By the freeness condi- 
tion (4), a non-empty periodic graph is necessarily infinite. A convenient represen- 
tation of a periodic graph in terms of a finite amount of data has been developed 
in a more general context in [12, 13], and probably independently in [3]. Due to 
the diversity of the literature spanning various fields of science, no universal termi- 
nology has been established. Here we partly try to follow the terminology of graph 
theory [14]. While this subsection contains standard material, we try to offer a 
slightly different point of view emphasizing the analogy to covering spaces [15, ch. 
1.3]. ' 

Given a group action on some mathematical object, it is natural to consider 
the quotient object with respect to the group action. For a periodic graph T, this 
means to identify two vertices (or edges) if they can be translated into each other 
by a group element x £ Z d ; in other words, if the two vertices (edges) lie in the 

same Z d -orbit. The resulting collection of vertex orbits V r /%d d = Vr/Z rf and the 

collection of edge orbits -Er/z d = £r/Z d form a quotient graph r/Z d : by (3), the 
source and target functions sr and ir descend to well-defined maps 

s r/z d ; tr/z d '■ Er/z d — ► ^r/z d • 
By the finiteness assumption of definition 2, the quotient graph r/Z d is finite. 

Remark 4. In the particular case that T C R d is a translation-invariant Euclidean 
periodic graph, then one can construct r/Z d also by taking the vertices and edges 
in a unit cell of T. Besides the edges inside the unit cell, each edge of T which 
connects a vertex inside the unit cell to a vertex outside the unit cell defines an 
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additional edge of r/Z d by changing the target vertex to its translate inside the 
unit cell. 

By definition of T/Z d , there is a canonical projection map (f>r ■ T — > T /Z d 
which maps every vertex and every edge to its Z d -orbit. When T is clear from the 
context, we also simply write 4> for <f>r- 

The map <fi enjoys the nice property that an edge (or a path) in r/Z d can 
be uniquely lifted to an edge (a path) in T, given that a starting vertex has been 
specified: 

Lemma 5. (a) For every e £ £r/z d an d every v* £ </> _1 (s(e)) ; there is a 
unique e* £ Er with </>(e*) = e and s(e*) = v* . 
(b) For every path p = e± . . . e n in T ' j1 d and every v* £ <p~ 1 (s(ei)), there is a 
unique path p* = e\ . . . e* in T with <j>(p*) = p and s(e*) = v* . 

Proof. (a) Let e £ Er be some edge with (f>(e) = e. Then there is a 

unique x £ Z d such that s(e) + x = v*. Hence, e* = f e + x has the 
desired properties. For uniqueness, suppose that e' £ Er would also satisfy 
<j>(e!) = e and s(e') = v*. By definition of (f>, the relation 0(e*) = e = <fi(e') 
means that e* and e' lie in the same Z d -orbit, so that there exists y £ Z d 
with e' = e* + y. But then, v* — s(e') = s(e*) + y = v* + y, which implies 
y = by (4). Hence e 1 = e*. 
(b) This follows from a successive application of part (a) to each edge in the 
path. 

□ 

In the language of graph theory, we have found that the projection tf> : T — > 
r/Z d is a covering of graphs [14, ch. 2], [1, ch. 17]. This is completely analogous to 
the notion of covering space in topology [15, ch. 1.3]. 

Unfortunately, knowing the quotient graph T /1 d alone is not enough to recon- 
struct r. For example, there are many Euclidean periodic graphs T C R d which 
contain only a single vertex per unit cell, so that iVp/z^l = 1- I n this case, all 
edges of r/Z d are loops. However, knowing \E T /%d \ as the number of these loops 
is certainly not enough to recover T: for example, it is unclear whether a loop of 
r/Z d comes from an orbit of loops in T, or whether it represents a class of edges 
connecting different vertices in T. 

The additional piece of data needed in order to recover T turns out to consist 
of edge weights on r/Z d with values in Z d , which specify, intuitively speaking, the 
translation required in going from Sp(e) to t r (e). These edge weights are known 
under various names — e.g. voltage assignments [12, 13], labels [3, 17], or simply 
weights [4] . We will prefer the term displacements, which we deem most appropriate 
given the geometric intuition. The formalism described in the following is closely 
analogous to gauge theory [18]. 

Defining displacements requires that one has chosen a vertex representative 
i(v) £ Vr for every orbit v £ Vp/z d - in other words, we fix a map I : Vp/z d Vr 
which is assumed to be a section of <f> : Vr — > ijfi . For example for a Euclidean 
periodic graph T C R d , one possibility is to define i by choosing a unit cell and 
mapping every orbit v £ Vp ^ d to its representative in the unit cell. 
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Every v G Vpma represents a whole 1 d worth of vertices of T, namely cf> 1 (v). 
Given l, we can choose a concrete identification of this </> -1 (v) with Z d by defining 

a v : Z d ^> <^ _1 (w) , x i — ^ f,(u) + x . 

This map is compatible with translations in the sense that it satisfies the identity 
a v (x + y) = a v (x) + y for all x, y G Z d . 

Now for an edge orbit e G -Ep /z d i lemma 5 provides an isomorphism (s(e)) J= > 
(i(e)), which is also compatible with translations. In total, we obtain an iso- 
morphism 

7e : Z d ^ r r 1 (s(e))^-r 1 We))^r-Z rf (5) 

which is again compatible with translations, r y e (x + y) = le(x) + y. Therefore, 

7 e (x) = 7 e (0) + a;. We now define the displacement along e to be 8(e) = f 7 e (0). The 
equation 7 e (0) = 8(e) expresses the intuition that 8(e) is the physical displacement 
required in going from Sr(e') to <r(e') for any e' G _1 (e). 

It is not difficult to see that the quotient graph T/Z d together with the dis- 
placement function 8 : E r / Z d — > Z d is sufficient to recover T. This works as in the 
following definition, which can also be regarded as a general scheme for constructing 
periodic graphs: 

Definition 6 ([13, 3, 4]). A displacement graph (G,8) is a finite graph G together 
with edge weights 8 : Eq — > Z d (the displacements ). Associated to (G,8) is a peri- 
odic graph G given by 

V s = {(v,x) | v G Vgi x G Z d } 
E d ={(e,x) | e G Eg, x G Z d } 

s c(( e > x )) = ( s o(e),x) , t r ((e, x)) = (t G (e) + 8(e), x) 

Remark 7. Intuitively, G is constructed from G as follows: we start with the 
lattice Z d and place at each point a copy of the vertex set Vq. Each edge e G Eq 
defines a Z d worth of edges between the copies of sc(e) and t G (e), where in adding 
these edges we have to apply a translation by 8(e) in the ambient Z d . 

One needs to keep in mind that determining a displacement graph from a 
periodic graph T requires choosing a representative l(v) G Vr for each Z d -orbit v G 
V r ijd . How much does the displacement graph depend on the choice of l1 Suppose 
we are given two different choices i, l' : Vr/z d — * ^r- Then for any v G Vr/z^i there 
is a unique g(v) G Z d such that 

l'(v) = i(v) + g(v) Vw G V r/Z d . 

By (5), the corresponding displacements 8 and 8' therefore differ by 

8'(e) = 8(e) + g(s(e)) - g(t(e)) Ve G E T/Zd . 

Example 8. Consider the periodic graph illustrated in figure 1(a). Any one of the 
elementary parallelograms formed by the dashed lines can be taken as a unit cell. 
Choosing the two vertices inside such a unit cell defines i in terms of a represen- 
tative of the "•" vertices which form a Z d -orbit, and a representative of the "■" 
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vertices which form another Z d -orbit. The associated displacement graph is shown 
in figure 1(b). One obtains the displacements of e.g. the edges going from • to ■ by- 
noting that there are three Z d -orbits of such edges in figure 1(a): one which stays 
inside the unit cell, corresponding to the displacement (0,0); one whose target is 
one cell away in the positive y-direction, having displacement (0, 1); and one whose 
target is one cell away in the negative x-direction with displacement (—1,0). 

If one would start from the displacement graph 1(b), one would probably draw 
its associated periodic graph as in figure 1(c), which is a different Euclidean em- 
bedding of the same periodic graph as in figure 1(a). 



3. Velocity Polytopes 

From now on, we assume T to be a periodic graph equipped with a fixed choice 
of orbit representatives b : Vr/z d —^Vr- If p = &\ ■ ■ ■ e n is a path in T/Z d , then by 
abuse of notation we define its displacement to be given by 

(6) 

i 

which is nicely coherent with the lifting properties of lemma 5: if each 6j lifts to 
an edge which intuitively translates by <5(ej), then the path e± . . . e n should lift to 
a path which intuitively tranlates by S(ei). 

3.1. Velocity. We now formalize the concepts introduced in the introduction. 
For technical reasons, we formally define a trajectory as a sequence of edges rather 
than vertices: 

Definition 9. A trajectory in T is a sequence (f n )neN of edges f n £ Er such that 
s(/n+i) = t(f n ) for all neN. 

Intuitively, a trajectory is nothing but an infinite path in T. Due to lemma 5, 
up to an overall translation a trajectory in T is uniquely specified by the sequence 
of edges (f>(f n ) S E r / Z d, which are the images under the projection 4> : T — > T/Z d . 
In the following, we will abuse notation by also writing /„ for <j){f n ). 

By the definition (6), the displacement traversed by the trajectory / between 
n = rt i and n = n%, i.e. along the path f ni . . . fn 2 -i, is given by 

Tl2— 1 

£ s (fk) ■ 

k—ni 

Since this displacement gets traversed in n,2 — n\ timesteps, it makes sense to define 
the velocity in that time interval to be given by the difference quotient 

n 2 - ni 

The trajectory / has a well-defined velocity if the limit 

= lim nu e Rd (8) 

n~ > oo Ji 

exists. In this case, the difference quotient (7) also converges to Uf for n% — > oo for 
any m G N. 

In the following, || • || will be a fixed but arbitrary norm on M. d . 
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(c) A different embedding of (a). (d) The nonzero basic velocities 

(arrows) and the velocity polytope 
(hexagon) . 

Figure 1 . A periodic graph in two different embeddings (a) , (c) , 
its associated displacement graph (b) and its velocity polytope (d). 
For more detail, see examples 8 and 13. 



Lemma 10. The velocity Uf of the trajectory f does not depend on the particular 
choice of representatives i : Vryz d — ^ used for constructing the displacement 
function S, but only on T and f themselves. 
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Proof. Let t, i! : Vp/z d —> Vr be two choices of orbit representatives. It has 
been noted in section 2.3 that their displacement functions satisfy 

S'(e) = 5(e) + g(s(e))-g(t(e)) 

for some appropriate function g : Vr/z d ~^ Then the displacements associated 
to the path f± . . . f n differ by 

n 

5'{h ■ ■ • fn) = S(h .../„) + E [9Wk)) - 9Wk))] ■ (9) 

k=l 

Due to s(/fe + i) = t(fk), the sum is telescoping, so that 

S'(h .../„) = S(h . . . f n ) + g{s(h)) - g(t(f n )) . 
Writing C = max^gy ||<7(v)||, we conclude that 



2C 
< — 

n 

from which the assertion immediately follows by taking the limit n — > oo. □ 

Note that the velocity of a finite path as in (7) is in general not well-defined in 
the sense of the lemma, but does depend on the choice of l. 

Remark 11. When r C R d is a Euclidean periodic graph, then we would like this 
notion of velocity to correspond to the usual one familiar from classical mechanics. 
As already mentioned in remark 3, we can always apply an affine transformation 
to r such that the unit cell becomes the ordinary unit cube [0, l] d . As a matter of 
bookkeeping, this also changes all velocity vectors by the same affine transforma- 
tion. 

We claim that if the unit cell of T is the unit cube, then the definition (8) 
gives precisely the usual concept of velocity. To see this, let us choose the orbit 
representatives i to be those in the unit cell, and write M for the maximal distance 
between any two vertices in the unit cell. Then, the actual distance vector traversed 
along the path /i . . . /„ will differ from the displacement J^k ^(/fe) by at most 2M. 
In the limit asn-> oo, this is negligible, since all distances get divided by the total 
elapsed time n. This proves the claim. 

As a trivial example, a constant trajectory has a velocity of 0. Similarly for 
any trajectory which stays in a bounded region in T C M. d . 

More non-trivial examples of velocities for not necessarily Euclidean T are as 
follows: for a cycle c = e\ . . . e n in r/Z d , we define the basic velocity associated to 
c to be given by 

Uc= s M = ELi 

n n 

This coincides with (7). Lemma 5 guarantees that c lifts to a unique path in L 
upon choosing an arbitrary vertex in </> -1 (s(ei)) as starting point. The trajectory 
defined by lifting a periodic traversal of c from r/Z d to T has the basic velocity u c 
as its velocity. 

One can regard the set of basic velocities as an invariant of the periodic graph: 

Lemma 12. For a given cycle c, the basic velocity u c does not depend on the 
particular choice of i : Vp i%& —¥ Vr used in constructing the displacement function 
S. 
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Proof. Applying equation (9) in this case, one finds that all terms in the sum 
cancel each other, so that 8'(c) = 8(c). □ 

Since by lemma l.(b) there are only a finite number of cycles in r/Z d , there is 
only a finite number of basic velocities for fixed T. 

Example 13. We go back to the periodic graph illustrated in figure 1. There are 
9 cycles in the displacement graph of figure 1(b) which have "•" as their starting 
vertex, all of length 2; the other 9 cycles with "■" as their starting vertex have the 
same basic velocities, so it is sufficient to consider the former. The basic velocities 
are 



(0,0) + (0,0) (0,0) + (0,-1) (0,0) + (1,0) 



2 






2 




2 


(0,1) + 
2 


(0,0) 


(0,1) 


+ (o,-i) 

2 


(0,1) 


+ (1,0) 
2 


(-1,0) + 


(0,0) 


(-1,0) 


+ (o,-i) 


(-1,0) 


+ (1,0) 



The nonzero ones are depicted in figure 1(d). 

Now that we have seen some examples of velocities, a natural question to ask 
is the following: 

Question 14. Given a periodic graph T, what is the set of its velocities? 

3.2. Main theorem. We now proceed to answer question 14 and give an 
explicit description of the set of velocities of T as a subset of M. d . We still take the 
periodic graph T with its associated displacement graph (r/Z d , 8) fixed. As before, 
I II will be a fixed but arbitrary norm on R d . The constant 

C = max 1 1 8(e) | 

will be of some use. |V| will always stand for |Vp/ Z d|. 

Lemma 15. For every path p in T/Z d , there are cycles Ci,...,Cfc in T/Z d such 
that 

fc 

<C|V| and < |p| ~^2\ci\ < \V\ . 



S(p) 



fc 



i=l 



8{ Ci ) 



Proof. For \p\ < \V\, there is nothing to prove since one can just take fc = 0, 
i.e. the sum over cycles to be empty. For \p\ > \V\, we use induction on |p|. By 
lemma l.(a), the path p can be written in the form 

p = p cpi 

where c is a cycle, so that the paths po and pi can be composed to p' = popi- Since 
|c| > 1, we conclude \p'\ = \p$\ + \pi\ < \p\, so that an application of the induction 
assumption to p' gives cycles eg, . . . , with 



fc 

8(p')-J2s(c i ) 

i=2 



< c\v\ 



id < \p\ - 



fc 

E 

i=2 



<\v\ 
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The conclusion follows by setting ci = c and observing 5(p) = S(ci) + 8(p') and 

IpI = |ci| + \ P '\. ' □ 

Theorem 16 (Main theorem). Let T be such that T /If 1 is strongly connected. Then 
the set of velocities of T coincides with the convex hull 

Pp = conv {y,f | Uf basic velocity in r/Z d } . (10) 

In particular, Pr is a rational poly tope in M. d , the velocity polytope ofT. 

We now give an outline of the proof before diving into the details. The idea 
is as follows: if Ci and C2 are cycles in T/Z d with the same initial vertex, then the 
closed path c\c^ has a velocity which is a convex combination of the basic velocities 
associated to c\ and C2. An analogous statement holds for longer combinations of 
cycles. 

So in order to show that the velocity of a trajectory is always a convex com- 
bination of basic velocities, one can apply lemma 15 in order to decompose the 
trajectory into cycles, noting that the right-hand side becomes irrelevant in the 
limit. 

Conversely, for every convex combination of basic velocities one needs to con- 
struct a trajectory which has this velocity. By choosing the number of times that 
each cycle appears in a closed path, one can adjust the coefficients of the convex 
combination which corresponds to the velocity of (the lift of) that closed path. 
Therefore, one can try to combine the cycles such that they appear in the tra- 
jectory with the appropriate frequencies, while also inserting some auxiliary paths 
which connect between cycles with different starting vertices. 

PROOF of Theorem 16. Working with the displacement graph (T/Z d ,S) in- 
stead of the periodic graph T, we begin by showing that any velocity lies in the 
convex hull of the basic velocities. We first claim that for any trajectory / and any 
b£N, the quotient 

n 

has the property that there is a vector u in the convex hull of the basic velocity 
vectors such that 

, 2\V\C , , 

\\w n -u \\< 1 ' 11 
n 

To this end, we first approximate the path /i . . . f n as in lemma 15 by cycles 
ci , . . . , Cfe , so that 



i=i j=i 



<C\V\ and < n - E l c il ^ \ V \ • ( 12 ) 

3=1 



The first inequality of these two can be rewritten as 



E 

3 



6(c 3 



II 



< 



C\V\ 



(13) 



Since the fractions S(cj)/\cj\ are basic velocities, the expression 
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is a convex combination of basic velocities; by the second inequality of (12), this 

— ■ appearing in (13). More precisely, 



approximates well the term ■ 



u- 1 - 



(13) C\V\ 
< — ! — - 



1 - 



Since ||u|| < C and n — ^ z |c/| < we can also bound the second term on the 
right-hand side by C\V\/n, which proves the claim (11). 

It follows from (11) that when the trajectory has a well-defined velocity lim„ w ni 
then the distance from this velocity to the convex hull of basic velocities is smaller 
than 2\V\C/n for any n. Since that convex hull is a polytope and therefore closed, 
the limit lim„ w n is itself in the convex hull of basic velocities. 

Conversely, it has to be shown that any convex combination of basic velocities 



(14) 



for weights Ai,...,A r > with J^Aj = 1 and cycles Ci can be realized by a 
trajectory. In order to construct such a trajectory, let us choose natural numbers 
ctik giving rational approximations to the numbers Aj/|cj| as 



Oiik 



clef 



k .2± 

M. 



so that 



Oik 

k 



< - Vi,keN. 
k 



(15) 



Let us also choose paths pi, . . . ,p r mT/Z d such that \pi\ < \V\ and pt connects the 
ending vertex of the cycle c, to the starting vertex of the cycle c,+i (with c r +i = c%, 
such that p r connects c r to c%). Such a choice of paths exists due to the assumption 
of strong connectivity. Then the building blocks of the trajectory are going to be 
the paths 

q k = f Ci . . . Cipi C 2 . . . C 2 P2 ...C r ...CrPv (16) 



ct\k times 02^ times 



a r k times 



Each q k is a closed path in T/Z d in the sense that its starting vertex coincides with 
its ending vertex. Its length can be estimated as 



\q k \ = ^2a ik \ci\ + ^2 \pi 

i—l i—1 i—1 



E 



k ■ 



A, 



c t \ + 0(l)) =k + 0(l) 



(17) 



where O(l) refers to a term which does not depend on k. This implies the similar 
estimate ||5(<Zk)|| < O(k). 

The trajectory / is defined to be (the lift to T) of the infinite path 

/ : glg2g2g393g3 ■ ■ • gfc ■ • ■ ■ ■ ■ ■ (18) 

k times 

It needs to be shown that this trajectory has a velocity which equals (14). We 
consider this trajectory up to a timestep n, i.e. the path f\...f n . We choose k 
(as a function of n) such that the path f\. . . f n has already concluded the whole 
gfc-segment from (18), but not yet the whole -segment. Then due to (17), 
the qk+i -segment will contribute to the time n and the displacement S by at most 
0(k 2 ). Hence the total displacement up to n timesteps is given by 

n k . a k r 
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Using condition (15) in the form ct!j m — < 0(1) evaluates this to 

i=l m — 1 1 VII / i— i 

On the other hand, the number n of edges traversed, which equals the time elapsed, 
can be evaluated very similarly as 

k k r 

n= E m ' Isml + o(fc 2 ) = E E m ( aim l c *l + + 0( fc2 ) 

m— 1 m— 1 i—1 

= E E w ( A * m + + °^ 2 ) = y E A * + °( fc2 ) = y + °( fc2 ) • 

m— 1 i—1 i=l 

Hence the velocity of the trajectory is given by 

Bm SkM) _ llm fgU^+og) _ £ *w 

n->oo n fc^oo ^.+0(^2) ^ \a\ 

as desired. □ 

Example 17. For the periodic graph of figure 1(a) and 1(c), the velocity polytope 
is the hexagon in figure 1(d). 

We now consider what happens when the connectedness assumption of theo- 
rem 16 is dropped. 

Proposition 18. Let T be any periodic graph. Let TifZ d , . . . , T c /1 d be the strongly 
connected components ofT/1 d which have preimages r 1; . . . ,r c C V. In this case, 
the set of velocities of T is the union of polytopes 

P r = P ri U...UPr c . (19) 

Proof. The I\ are defined as the preimages under cj) : T -> T/Z d of the 
strongly connected components of F/Z d . Then for every trajectory (f n )neN, there 
is no G N such that all /„ for n > no lie in the same Fi . The velocity of / therefore 
lies in the corresponding velocity polytope Pr r This proves the "C" inclusion 
of (19). The "3" inclusion is clear since any trajectory in Ti is also a trajectory in 

r. □ 

4. The large-scale geometry of a periodic graph 

What we mean by the large-scale geometry of a strongly connected periodic 
graph r is the following. There is a natural notion of distance between two vertices 
defined to be the length of the shortest path connecting the two vertices. This 
defines a metric d(-, •) on T invariant under the action of % d . Given any x G Z d and 
any vertex ti £ l^r, we can now define the T-norm ||x||r to be given by 

||x|| r d = f hm + (20 ) 

The existence of the limit is guaranteed by the triangle inequality and Fekete's 
lemma. Taking the limit instead of defining | \x | | r to be d(v,v + x) itself is necessary 
in order guarantee that ||x||r does not depend on v. This well-definedness of (20) 
easily follows from arguments very similar to those made in the proof of lemma 10. 
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Moreover, it is simple to show that ||mx||r = m IM|r for any m £ N. The norm 
||x||r can be interpreted as follows: for any ngT, one needs to traverse n ■ \ \x\\r + 
0(1) edges in order to get from v to v + xn. 

From the triangle inequality \\x + y\\r < \\%\\r + IMIr & n d ||mx||r = m||x||r 
for m £ N, one deduces that || • ||r extends to a unique norm on M. d . 

Proposition 19. Let T be strongly connected. Then Pr is the unit ball of \ \ ■ \ \-p- 

Proof. We first show that Pr is contained in the unit ball. To this end, it is 
enough to prove ||w c ||r < 1 for & basic velocity u c associated to a cycle c = e± . . . 
in r/Z d . Choosing any v £ </> _1 (s(ei)), the cycle c lifts to a path in V from v to 
v + 5(c). Since this path has length k, we get 



Klli 



5(c) 



< \ lim ^ 



+ n5(c)) I u 

< -rd(v, v + 5(c)) 

n k 



There is a path from v to v + 5(c) of length fc, which implies d(v, v + 5(c)) < k. 
This results in the desired inequality ||w c ||r < 1- 

Conversely, it has to be shown that ||x||r < k for x € Z d and k € N implies 
that f € Pr- Fix e > 0. By assumption, there is v £ T and n £ N such that 



e < 



d(v, v + nx) 



<k + e. 



This means that there exists a path p from v to v + nx of length between nk — ne 
and nk + ne. Since the final vertex of p is a translate of its starting vertex, p 
can be concatenated with its own translates in order to form a trajectory which 
periodically traverses translates of p. Since 5(p) — nx, the velocity of this trajectory 
is given by nx/\p\ 7 so that 

x kn _ 
G -Pr • 



In terms of an arbitrary norm 
bounded by 

kn 



k \p\ 

| on R d , the distance of f to Pr can therefore be 



x 
k 



\P\ 



k 



1 



|p| — nfc 



< 



Fl 



Since this vanishes as £ — J> 0, we conclude that f G Pr from closedness of Pp 



□ 



5. Properties of velocity polytopes 

We now study some basic properties of velocity polytopes. 

In section 3, we have assigned a velocity polytope to each (reasonably con- 
nected) periodic graph. By remark 11, for a Euclidean periodic graph T C R d , this 
polytope describes the set of achievable velocities of a classical particle hopping 
along the edges of T. In order to get interesting models for physics, the velocity 
polytope should be not too different from a Euclidean ball, thereby making the set 
of achievable velocities close to isotropic. This raises the question, which rational 
polytopes can be realized as velocity polytopes of periodic graphs? This is not 
difficult to answer: 

Proposition 20. Every rational polytope arises as the velocity polytope of an ap- 
propriate periodic graph T. 
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PROOF. Let P C M. d be a non-empty polytope with rational vertices w\ , . . . , w m 
Q d . Let 7 be the least common multiple of the denominators of the Wi, so that 

G 7L d for all i. Then P can be realized as a velocity polytope as follows: let 
us construct a displacement graph on vertices u\, . . . ,u 7 such that there is a single 
edge each with displacement <5 = from the vertex Uj to the vertex Uj + \ for each 
j = 1, . . . , 7 — 1, and additional edges e\, . . . , e m from u 7 to u\ with displacements 

5(ei) =' ^ . 

This defines a strongly connected displacement graph. Its basic velocities are pre- 
cisely all the — . Therefore, the associated periodic graph has P as its velocity 
polytope. □ 

In this paper, we have been considering the general case of directed graphs; 
since every undirected graph can be made into a directed graph by replacing an 
undirected edge by a pair of oppositely oriented directed edges, the theory also 
applies to undirected graphs. In an undirected graph, every path can be reversed, 
which reverses the sign of its velocity. Therefore, it should not be surprising that 
the following holds: 

Proposition 21. When T is undirected, then Pr is symmetric around the origin. 

Proof. The edges of T come in parallel pairs of opposite orientation. There- 
fore, the edges of r/Z d also come in parallel pairs with opposite orientation and 
displacement of the opposite sign. Hence for every cycle c in L/Z d , there is a cycle 
c' which corresponds to traversing c backwards by using all the "partner" edges. 
Therefore, if u c is a basic velocity, then so is — u c . Now the statement follows from 
theorem 16 and proposition 18. □ 

Proposition 22 (Connectedness). IfT itself is strongly connected, then Pr C M. d 
is full- dimensional and G Pr is an interior point. 

Proof. This is best proven without appealing to theorem 16. Let ei, . . . , ed G 
1i d be the standard unit vectors, and v G T a fixed starting vertex. Then by the 
assumption of strong connectedness of T, there is a path pf in L from v to v + e,. 
In L/Z d , this is a closed path of displacement ti and velocity j^jrr. We can use 

translates of pf to connect v + nti to v + (n + l)ej for any n G N. Sequentially 
traversing these translates of pf defines a trajectory with velocity j^t . Similarly, 

there exists a path p~ from v to v — ei, which gives rise to a trajectory with velocity 
— t^t. The convex hull of these 2d velocities is a subset of Pr. By construction, 

this subset is full-dimensional and includes the origin as an interior point, and so 
the same also holds for Pr- □ 

The converse is not true: for example, in dimension d = 1 we can take Y fL d 
to consist of a single vertex with a loop e + of displacement S(e + ) = 2 and a 
loop e~ of displacement <5(e_) = —2, which makes T/Z d strongly connected and 
gives Pr = [—2,2], although the associated periodic graph L has two connected 
components. 

So far, we have been talking about the velocity polytope of a single periodic 
graph P. But given two periodic graphs T and V , how do their velocity polytopes 
relate? In order to find some relation between Pr and Pr', one needs to assume 



1(5 



TOBIAS FRITZ 



a relation between F and F' . One such notion is that of a morphism h :T —¥ V 
between periodic graphs V and F' of the same dimension d, by which we mean maps 

hy : Vr — ^ Vr' i h*E '• -Er — > E-p' , 

which are compatible with the graph structures, 

sr> o He — hy ° sr , *r' ° hE — hy ° tr , 

and the Z d -action, 

hv(v + x) — h v (v) + x , h E (e + x) = h E {e) + x Vv e V r , e eE T ,xeZ d . (21) 

A morphism of periodic graphs induces a simple relationship between the velocity 
polytopes: 

Proposition 23 (Functoriality) . Let h :F —±F' be a morphism of periodic graphs. 
Then 

Pr C Pr' • 

Proof. We consider how h operates on the quotient graphs F/Z d and F' jZ d . 
Since by (21) the assignment hy maps Z d -orbits to Z d -orbits, we get induced maps 
hy ■ V r /id — >■ V T : and h E : E r / Z d — ^ E T , / Z d . Moreover, the orbit representatives 
i : Vp/z d ~* Vr and i' ■ Vpi /jd Vr' can be chosen such that hy maps represen- 
tatives to representatives in the sense that hy o i = / o hy. In total, there is a 
commutative diagram 

r/z d -^r/z d 



such that the vertical compositions are identities. Then it follows that the induced 
map h E ■ E r /%d — > E r , /%d is compatible with the displacements, 5'(h(e)) = <5(e) for 
all edges ee£ r nt . 

To prove the assertion, we now show that every velocity of T, associated to a 
trajectory (f n )neN-> 1S also a velocity in F'. But this follows from 

lim ELiW = £LiW/fc)) ( 

n^oo n n 
so that (h(f n )) n£N is a trajectory in F' with the same velocity. □ 

This result can be applied for example as follows: 

Corollary 24. If Pr 2 -Pr'; then there is no morphism h : F — > F' . 

5.1. A no-go theorem for digital physics. We end the paper by concluding 
from theorem 16 and proposition 18 a no-go theorem for the construction of models 
in "digital physics" [10]. 

Theorem 25. There is no periodic graph F for which the set of macroscopic ve- 
locities achievable by a classical point particle hopping along the edges is isotropic. 
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Proof. In order to be isotropic, the set of velocities would have to be a Eu- 
clidean ball. However, a Euclidean ball is not a polytope or a union of a finite 
number of polytopes. □ 

References 

[1] Norman Biggs. Algebraic graph theory. Cambridge Mathematical Library. Cambridge Uni- 
versity Press, Cambridge, second edition, 1993. 

[2] Vladislav A. Blatov and Davide M. Proserpio. Periodic- Graph Approaches in Crystal Struc- 
ture Prediction, pages 1-28. Wiley- VCH, 2010. 

[3] S. J. Chung, Th. Hahn, and W. E. Klee. Nomenclature and generation of three-periodic nets: 
the vector method. Acta Crystallographica Section A, 40(l):42-50, 1984. 

[4] Edith Cohen and Nimrod Megiddo. Recognizing properties of periodic graphs. In Applied 
Geometry and Discrete Mathematics - The Victor Klee Festschrift. AMS, 1991. 

[5] L. Collatz. Spektren periodischer Graphen. Resultate Math., l(l):42-53, 1978. 

[6] Giacomo Mauro D'Ariano. The quantum field as a quantum computer, 2010. arXiv:1012.0756. 

[7] Giacomo Mauro D'Ariano and Alessandro Tosini. Emergence of space-time from topologically 
homogeneous causal networks, 2011. arXiv:1109.0118. 

[8] Olaf Delgado-Friedrichs and Michael O'Keeffe. Crystal nets as graphs: Terminology and 
definitions. Journal of Solid State Chemistry, 178(8):2480 - 2485, 2005. 

[9] Jean-Guillaume Eon. Euclidian embeddings of periodic nets: definition of a topologically 
induced complete set of geometric descriptors for crystal structures. Acta Crystallographica 
Section A, 67(l):68-86, Jan 2011. 
[10] Edward Fredkin. An informational process based on reversible universal cellular automata. 

Physica D: Nonlinear Phenomena, 45(l-3):254 - 270, 1990. 
[11] Chris Godsil. Periodic graphs. Electron. J. Combin., 18(l):Paper 23, 15, 2011. 
[12] Jonathan L. Gross. Voltage graphs. Discrete Math., 9:239-246, 1974. 

[13] Jonathan L. Gross and Thomas W. Tucker. Generating all graph coverings by permutation 

voltage assignments. Discrete Math., 18(3):273-283, 1977. 
[14] Jonathan L. Gross and Thomas W. Tucker. Topological graph theory. Dover Publications Inc., 

Mineola, NY, 2001. Reprint of the 1987 original. 
[15] Allen Hatcher. Algebraic topology. Cambridge University Press, Cambridge, 2002. 
[16] Franz Hofting and Egon Wanke. Minimum cost paths in periodic graphs. SIAM J. Comput., 

24(5):1051-1067, 1995. 

[17] W E Klee. Crystallographic nets and their quotient graphs. Crystal Research and Technology, 
39(11) :959-968, 2004. 

[18] George Svetlichny. Preparation for gauge theory, 1999. arXiv:math-ph/9902027. 

[19] Gcorg Thimm. Crystal topologies - the achievable and inevitable symmetries. Acta Crystal- 
lographica Section A, 65(3):213-226, May 2009. 

[20] Thomas Zaslavsky. A mathematical bibliography of signed and gain graphs and allied areas. 
Electron. J. Combin., 5:Dynamic Surveys 8, 124 pp. (electronic), 1998. 

ICFO - Institut de Ciences Fotoniques, Mediterranean Technology Park, 08860 
Castelldefels (Barcelona), Spain 

E-mail address: tobias.fritz@icfo.es 



